import sympy as sp
x, t = sp.symbols('x t')
H = x*sp.integrate(sp.sin(t**2), (t, 0, x)) - sp.integrate(t*sp.sin(t**2), (t, 0, x))
H_double_prime = sp.diff(H, x, 2)
print("H''(x) =", sp.simplify(H_double_prime))